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Abstract 

Light-cone QCD sum rules are employed to compute the strong coupling 
constants: gB*B*n, gBiBon an d gB^B*-*, where (B,B*) and (Bq,Bi) are nega¬ 
tive and positive parity (0~, 1 _ ) and (0 + , 1 + ) qQ doublets. The couplings are 
calculated both for finite values of the heavy quark mass and in the infinite 
heavy quark mass limit, deriving sum rules for tbq —> oo. 
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I. INTRODUCTION 


The recognition of the approximate symmetries of QCD, explicitely or spontaneously 
broken, represents a basic step for the understanding and the description of the strong 
interactions in processes involving hadrons. A remarkable example, in the sector of systems 
containing one heavy quark, is the SU(2Nf) heavy flavour-spin symmetry holding in the 
infinite heavy quark mass limit: rriQ f —> oo. In this limit, the masses of the Nf heavy 
flavoured quarks are irrelevant in the interactions with the hadronic light degrees of freedom, 
since the heavy quarks only act as static colour sources. Moreover, the decoupling of the 
gluon from the spin sq of the heavy quark allows to relate the properties of the states 
belonging to the doublets obtained combining sq with the spin s[ of the light degrees of 
freedom. The SU(2Nf ) spin-flavour symmetry is explicitely broken by the finite values 
of the heavy quark masses, and the symmetry-breaking effects can be described using an 
expansion in the inverse masses of the heavy quarks |I],[|]. 

At the opposite energy scale, for vanishing masses of the up, down and strange quarks, 
QCD is invariant under chiral SU(3)l x SU(3)r transformations; this symmetry is sponta¬ 
neously broken, the Goldstone bosons being represented by the octet of light pseudoscalar 
mesons [[|. 

Both the heavy quark spin-flavour and the chiral symmetries can be realized writing 
an effective QCD lagrangian, expressed in terms of hadronic fields 0. The interactions of 
the heavy sf — \ and | + mesons with the octet of light pseudoscalar 0“ mesons can be 
described by the lagrangian: 

£ = i Tr{H b v»D^ ba H a } + f -fTr{d» Ec^} + Tr{S b (i A‘D, lba - 5 ba A)S a } 

O 

+ i g Tr{H b j ll 'y 5 A ba H a } + i g' Tr{S b ^^ b A^ a S a } + [ih Tr{S b ^^ b A% a H a } + h.c] . (1.1) 


The fields H a in (|1 . 1|) describe the negative parity J p = (0 ,1 ) qQ doublets (with s[ — ^ ): 


H a = ( Ar L \P’ v Y - PM . 


( 1 . 2 ) 
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where the operators P/ M and P a respectively annihilate the 1“ (P*) and 0“ {B a ) mesons of 
four-velocity v (a = u,d,s is a light flavour index). In the infinite heavy quark mass limit 
such states are degenerate in mass, due to the vanishing of the chromomagnetic interaction 
of the heavy quark spin with the spin of the light degrees of freedom of the mesons. Hence, 
H a describes a doublet of states with the light degrees of freedom with zero orbital angular 
momentum with respect to the heavy quark. Analogously § , the fields S a describe a doublet 
of states having the light degrees of freedom with angular momentum (P-waves in 

the constituent quark model): 

S a = _ Do ] , (1.3) 


with Di, Dq annihilation operators for the states 1 + and 0 + , respectively. Notice that all the 
heavy held operators in the lagrangian (|1 . 1|) contain a factor y/mp and have dimension 3/2; 
the parameter A in (|1 . 1[ ) represents the mass splitting between the positive and negative 
parity states. 

The octet of light pseudoscalar mesons is included in the effective lagrangian (|1.1| ) using 

ij\4 r. 

the exponential representation: £ — e~T* and £ = £ , the matrix M. contains the n,K and 
7] fields: 


M = 


V 




7r 

K~ 


1] K + ^ 

K ° -\Hn) 


+ 


(1.4) 


with = 132 MeV. Finally, the operators D and A in (|1 . 1| ) are given by: 


Dfj,ba — fibadf! + V^ba ~ $bad^ + — (^8^ + ( 1 - 5 ) 

At>a = l(£d^-Zdy) ba . ( 1 . 6 ) 


The effective Lagrangian (|1 . 1| ) can be generalized to include the interactions of heavy 
sf = | + states, as well as the interaction with the octet of low-lying vector mesons [||]. 

The strong interactions between heavy H ai S a mesons and the light pseudoscalar mesons, 
as described by o, are determined by three couplings: g, h and g'. The spin symmetry 
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implies that the couplings to light mesons (from now on we refer to pions) do not depend on 
the particular member of a doublet; for example, the couplings gs*B*ir and Qb*b-k, defined 
respectively by the matrix elements 

< B*°(p,e 2 )n~(q)\B*-(p + q,e x ) >= -r 9 B*B^e prTOl pe p 1 e* 2 a p a q (i (1.7) 

J7V 

and 


2 m 


<B°(p) 7T (q)\B* (p + q,e) >= ^-^g B *Bne p q P 

J 7T 


( 1 . 8 ) 


(e is the B* polarization vector), coincide with the same coupling g in the limit mb —> oo. 
Analogously, the constant h describes the coupling to pions of whatever a member of the 
H a doublet and a member of the S a one, whereas the strong interaction of positive parity 
S a states with pions is governed by the coupling g'. 

Differences between the couplings of states belonging to the same doublet are of order 
1 jrriQ. Such differences play an important role in the heavy meson phenomenology; for exam¬ 
ple, the difference between the pion-vector-vector gs*B*iT and the pion-pseudoscalar-vector 
gB*Bn coupling constant appears in the 1/rriQ correcting terms to the effective lagrangian 
(O) 0, and must be taken into account in the calculation of the heavy meson hyperhne 
splitting 


In refs. |pT)| [T2[| the couplings g and h were computed using QCD sum rules and the short- 
distance expansion Jl3|], in the infinite heavy quark mass limit. The short-distance expansion 
was also adopted to calculate the couplings gs*Bin gB*B*-n and gsoB-K for finite masses of the 
heavy quarks |T^,|T0|,|TI|.|^.[T5[ ; comparing the results derived in the asymptotic regime and 
for finite masses, significant information was obtained on the size of 1/rriQ corrections. 

A different approach to compute the strong couplings was proposed in [JT6[] and adopted in 
TO, using sum rules based on the expansion near the light-cone [ITS EUJ Pi. The parameters 


1 Applications of the light-cone sum rule method to semileptonic heavy meson decays can be 
found in |2l| and in references therein. 
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9 b* B-m 9 Bo Bit were determined by this technique for finite masses of the charm and bottom 
quarks, and the couplings g and h were obtained extrapolating to rriq —> oo the results for 
finite m c and m 

In this paper we derive light-cone QCD sum rules in the limit rriq —> oo. This allows 
us to perform a systematic comparison between the results of light-cone and short-distance 
QCD sum rules, which can shed light on the critical parameters, the advantages and the 
drawbacks of both the approaches. As we shall see, in the infinite heavy quark mass limit 
a remarkable agreement between the two approaches is obtained, supporting our confidence 
on such methods for the calculation of parameters which are sensitive to the nonperturbative 
dynamics of QCD. 

The plan of the paper is as follows. In Sec. II we apply light-cone sum rules to calculate 
the coupling g and gB*B*ir ; the numerical analysis of this channel is presented in Sec. Ill, 
together with a comparison with the results obtained for gB*Bn■ In Sec. IV and Sec. V we 
analyze the couplings g' and h, respectively; in particular, we consider the process Bi —» B 0 n. 
In Sec. VI we draw our conclusions. 


II. LIGHT-CONE QCD SUM RULES FOR g B *B* 7 r AND g 


Let us first consider the finite heavy mass case, and compute the coupling gs*B*ir defined 


in eq. (|1.7|) . The correlator of two vector currents J M = q'y^b between the vacuum state and 
a pion external state with momentum q: 


Ffru(p,q) = i J dx <n (g)|T[d(x) 7 / , 6 (x), 6 ( 0 ) 7 I/ u( 0)]|0 > e ipx 
= ie^ a qP a q t3 F 


( 2 . 1 ) 


corresponds to the diagram depicted in Fig.l. Considering a pion on mass shell: q 2 = 

(we set m n = 0), the invariant function F depends on p 2 and (p + q) 2 . 

Following the same idea underlying the usual QCD sum rule method, a hadronic repre¬ 
sentation can be given to the correlator ( 0 ) in terms of the lowest lying resonance, the B* 
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pole, of higher resonances and a continuum of states. In a double dispersion relation for F 
in the variables p 2 and (p + q) 2 : 


F(p 2 , (p + q ) 2 ) = f dsds'- - 

J [s — 


p had (s , s') 


( 2 . 2 ) 


(s -p 2 )(s' - (p + q) 2 ) 

the spectral function p had gets contributions from B* and from higher states; therefore, it 
can be modeled as 


P had (s , s') = -r9B*B*Trm 2 B*fB*fi(s ~ m B *)S(s' - m|.) + p cont (s , s')0(s - s 0 )6>(s' - s' 0 ) (2.3) 

JIT 

where p cont (s , s') includes the contribution of the higher states and of the continuum above 
the effective thresholds so, s' 0 ; the B* leptonic constant f B * in ( p.3|) is defined as 


< 0\d r y lx b\B* o (p, e) >= m B *f B * 


'M • 


(2.4) 


We neglect in the dispersion relation (|2.2| ) possible subtraction terms, which play no role in 
the Borel transformed sum rules, as shown in the following. 

For large positive values of — p 2 and —(p+q) 2 the function F(p 2 , ( p+q ) 2 ) can be computed 
in QCD. This task is afforded expanding the T-product in eq.( [2.1| ) near the light-cone 
(. x 2 = 0) in terms of non local operators, whose matrix elements are defined by wave functions 


associated with operators of increasing twist. According to the notations in ref. |16|, the 
catalog of the relevant matrix elements, up to twist 4, is: 


< 7r(g)|d(a;)7^75M(0)|0 > = - if n q M / du e luqx ((p n (u) + x 2 g 1 (u) + 0(x 4 )) 

Jo 

2 i 

+ - —) [ du e iuqx g 2 {u) , (2.5) 

qx Jo 

— f 777 ^ /*! 

< 7r(g)|d(a;)i75M(0)|0 > = ———-— / du e luqx (p P {u) , (2.6) 

rn u + rrid Jo 

— f m 2 

< Tr(q)\d(x)a flu 'y 5 u(0)\0 > = i(q^x u - q v x M )——-—- / du e iuqx ip a (u) . (2.7) 

vyTTl u i TTld) J 0 

The wave function is associated with the leading twist 2 operator, g\ and g 2 correspond to 
twist 4 operators, and p> B and tp a to twist 3 ones. Due to the choice of the light-cone gauge 
x^A^ix) = 0, the path-ordered gauge factor P exp (ig s Jq dux^A^ux)) has been omitted. 
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Notice that the coefficient in front of the r.h.s. of eqs. ([2.6|) , ( [2.7| ) can be written in terms of 

f 7T 7Ji 2 2 

the light quark condensate < uu > using the PCAC relation: ———-— = —— < uu>. 

m u + rn d 

The pion matrix elements of quark-gluon operators can also be parameterized in terms 


of wave functions [If 


< 'K(q)\d(x)a a/3 'y 5 g s G^(ux)u(0)\0 >= 

'ifs-n [(.QfiQaQvfl QvQaQfif)) (Q[J,Ql39i'ot QvQ09fJ-a')\ / ^37r(c^i)c ^ ^ 


( 2 . 8 ) 


and 


< 7r(q)\d(x)'-/ ll 'y 5 g s G a i3(vx)u(0)\0 >= 
fnlqplgan - ^ £ ) - q a (g^ - ^ i )] / Va i (f ± (a i )e 

Q * X Q ’ X J 


D iqx{ai+va3) 


% 


+U——(q a X/3 - qpx a ) / T>ai<pn(ai)e 


l iqx(a 1 +vo3) 


q ■ X 


(2.9) 


< n(q)\d(x)j IJ ,g s G a i 3 (vx)u(0)\0 >= 


ifir[qp(. 


gan 


x a q^\ ( 


q ■ x 


T>a>i<p±(ai)e 


iqx(ai+vois) 


+ifn-^-(q a xp - qpx a ) I Va^^a^e 


iqx(a.\-\-voL 2 >) 


q • x 


( 2 . 10 ) 


The operator G ai g is the dual of G a p: G a/ g = \e a ps P G Sp ; Dcy is defined as Va>i = 
da\da2da^5{l — oq — a 2 — 07 ). The function 933 ^. is twist 3, while all the wave functions ap¬ 
pearing in eqs.( [2.9|) , ( |2.10| ) are twist 4. The wave functions </?(oy, /i) (// is the renormalization 
point) describe the distribution in longitudinal momenta inside the pion, the parameters .77 
(J2i x i — 1 ) representing the fractions of the longitudinal momentum carried by the quark 
and the antiquark. 

The wave function normalizations immediately follow from the definitions ( p.5| )-( pT0|) : 
/0 du ip n {u) = /g 1 du ip<j{u ) = 1, /g 1 du g±(u) = 5 2 /12, fd du G 2 (u ) = <5 2 /18 (with G 2 dehned 
as G 2 (u) = - Jo dv g 2 {y)) JVai<p ± (oii) = f Va^w (cy) = 0, J Dawson) = - J Va^on) = 
<5 2 /3, with the parameter 5 defined by the matrix element: < 7 r(g)|d 0 s (j a/i 7 “M|O >= id 2 f^q^. 

I 11 the limit g —> 00 the wave functions assume an asymptotic form which is approxi¬ 
mately symmetric in the variables ay and hence represents an almost equal distribution of the 
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pion momentum between its constituents. This asymptotic wave function can be calculated 
by exploiting the conformal invariance of QCD in the short distance region. On the other 
hand, to investigate the general properties of the wave functions and their deviation from 
the asymptotic form, non perturbative methods are required. Usually, the momenta of the 
wave functions are determined by QCD sum rules, and then the wave functions themselves 
are reconstructed [|16|,^(]]. In a calculation of the type presented here, however, the complete 
form of the pion wave functions is not required, since only the values in particular points, 
namely near the symmetric point uq = 1/2, are needed. Therefore, the uncertainty in the 
results is only related to a finite number of input parameters, as it will be shown in the 
following. 


The expansion of the correlator ([2.1|) near the light-cone produces the following expression 
for the invariant function F: 


F QCD (p 2 , (p + q) 2 ) = f 

Jo 


du 




m 2 — (p + uq) 2 




du 


o [ml — (p + uq) 2 ] 2 L ml — (p + uq) 2 


1 + 


2 ml 


(gi(u) + G 2 {u)) 


2771b- 


Uml 


r-l 


du-, 


ip a {u) 


'Q{m u + m d )Jo [ml - (p + uq ) 2 ] 2 

, f 1 , /' yuK) -y±(Q!») 

+ U Jo J '{ml-lp + qia.+vas )] 2 } 2 

r 1 r i 


( 2 . 11 ) 


p-q 


4 U dv(v-l) da^{a 3 )^- t— 

Jo Jo [ml - [p + < 7(1 + (v - l)a 3 )pu 

f 1 /•! rl—a3 „ p . q 

4/tt dv da 3 / d«i 0 (ai, a 3 )-p —^ -f—i —7 -;-vPnT 

Jo Jo Jo {ml — [p + q{a \ + vas)] 2 }- 5 


+ 2 J v dv J Vai 


^||(«* 


/o J ‘ {m 2 b - [p + g(«i + ra 3 )] 2 } 2 
where 7p(af) = — / 0 “ 3 dt / 0 1_ * 0 («i, t)da\, 0 (ai,a 3 ) = — / 0 ai dt 0 (t, a 3 ), and (j) = (p± — <p± + 
<p\\ — (p\\. Taking into account the wave function normalizations we recover, in the soft pion 
limit q —■> 0 , the expressions obtained in H], 


The sum rule for gB*B*n follows from the approximate equality of the expressions ( |2.2| ) 
and ( |2.11| ). Moreover, invoking duality arguments, the contribution of the continuum in 
( p.2|) corresponds to the QCD contribution. This allows to isolate the pole contribution to 




















( j2.2|) . The equality between Q2.2|) and (|2.11|) can be improved performing a double Borel 
transform in the variables —p 2 and ~(p + q) 2 - Defining M 2 and M 2 as the Borel parameters 
associated to the channels p and p + q, respectively, and using the formula: 

(*-!)! 


ft M\ Bm 2 


2 [m 2 — (p + uq) 2 ] 1 M 2 M 2 

M 2 


(■ M 2 ) 2 ~ e , m 2 + q 2 (l-u 0 ),^ 

exp (- jp -) S{u - u o> 


( 2 . 12 ) 


. , , , 2 M(M£ 
with M 2 = ^-ttv and Uq = 


Ml + Ml 


Ml + Ml 


, we get: 


F QCD (Ml,Ml) = 2 e M 2 (p n (u 0 ) + 


m b mt 


M{M. 


3 (m u + m d ) 


<Ar(w 0 ) 


~ + ~m?) ^9i( u o) + G 2 (u 0 )] 

1_Q1 da 3 


ru o 


-b / dcxi 

J 0 JuQ—ai 0^3 


(ip ± - <f_ l) (cq, 1 - Oil - OL 3 , « 3 ) 




dan 


f 1 —ai 


da ? 


«0 — Cli 


/WQ — 




< / J ||( a l; 1 ~ Q!l — Q!3, CI3) 


(2.13) 


$(! - «o) 


den 


^(«3)l 


/l-n 0 




L Jo 


1 - Wo 

r° da 4 >(wo - «3,«3) r 1 da 0(w o ,a 3 )i 

CI3 do «3 


On the other hand the Borel transformed hadronic representation of the function F reads: 
2 m 2 B »fl»g B * B *TT 


F(M*,Mi) = 


fn Ml Ml 


T ^ + / dsds'p con \s, s')e "T "2 


(2.14) 


As shown in ( p,14|) , the Borel transformation exponentially suppresses, for small values of 
the parameters Mf, Ml, the contribution of higher states and of the continuum; moreover, 
possible subtraction terms in Q 2 . 2 |) , which depend only on p 2 or (p + q) 2 , are removed by 
the independent borelization in the two channels. This is a feature which renders light-cone 
sum rules quite appealing, as far as the calculation of the strong couplings is concerned. 
As a matter of fact, the determination of the strong couplings by, e.g., short-distance rules 
is usually performed in the soft pion limit q —> 0 , thus preventing the possibility of an 
independent borelization in the two B* channels. 

Considering the symmetry of the correlator ( |2.1| ) (Fig.l), it is natural to choose M( = 
Mf, and then w 0 = 1/2, together with s 0 = Sq- Such a choice corresponds to a quark and an 
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antiquark of the same momentum inside the pion. In this condition the subtraction of the 
continuum can be performed, in the leading twist term, with the substitution e~ m V M2 —> 
e ~ m l/ M 2 _ e -so/M 2 a rec jp e we follow considering that the numerical contribution of the 

higher twist terms is small. 

In this way we get the sum rule for gs*B*F- 


<■2 

J b*9b*b*k — 


f- 

J 7 T 


2 


e m* l (e m 7 — e m 7 )M 2 (p n (uo) 


+ e m 7 

ru o 


+ 


da 


^ ( “ o) “ 4 (' + lid [9l( “ o) + G2( “ o)1 

1- “i da 3 


lo 


/ UQ—ai «3 


{<P± ~ <P±) («i, 1 - «i - a 3 , a 3 ) 


f u ° rl—ai Uq _ q , 1 

2/ da i/ da 3 -,— Vdl («i, 1 - «i - «3, “ 3 ) 

JO Juo-a-L a 3 


(2.15) 


r$(l - '“ 0 ) 


fl 


den 


^(«3)l 


Mo 


' 1 —no 


C^q 


f u ° da 0 (m o - a 3 ,« 3 ) 

L Jo a 3 


/o a 3 


The numerical analysis of this sum rule is described in the next section. Here we want to 
show that eq. (|2.15|) has the right behaviour in the limit mj —> 00 jHj , and that an analytical 
expression can be derived in the asymptotic regime. 

In the limit mb —»■ 00 the B* meson mass can be related to the b quark mass by the relation 
ms* = mft+A + Ch—), where A represents the binding energy of the light degrees of freedom 
in the static b quark chromomagnetic held. Moreover, the B* lcptonic constant depends 
on rrib as: /#* = -y= (modulo logarithmic corrections), with F a low energy parameter 
remaining constant in the infinite limit. Rescaling the Borel parameter: M 2 = 2mbE , 
and the continuum threshold s 0 = m l + ^ m byo-, with E and y$ independent of the heavy 
mass, we recover from (|2.15|) a sum rule holding in the asymptotic limit: 


F 2 g = SF |(l ")EiPn(u„) - < “J 2 > <p„(u 0 ) - -tjfol(“ 0 ) +G 2 (“o)1 


(2.16) 


which is the analogous of the sum rule derived for g 0 by the short distance expansion. 
Notice that the sum rule ( j2.16| ) is not hampered by the problem of the subtraction of the 
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so-called ’’parasitic” contributions, which must be carefully treated in the short-distance 


calculation [|T 0 |; as a matter of fact, it results from two independent borclizations with 
parameters E\ and E 2l with u 0 = Ei+e 2 anc ^ ^ = Ei+e 2 • 

Finally, for mb —^ cxo, it is worth observing that the light-cone sum rule for Qb*b-k |16 
coincides with eq.(|2.16|), as required by the heavy spin symmetry. 


III. NUMERICAL ANALYSIS FOR g B *B*n AND g 


The numerical analysis of eqs. (|2 .1 5|) , (|2.16| ) requires a set of light-cone pion wave 
functions. Actually, if the choice u 0 = - is made, only the value of the wave functions in 
the middle point is needed. 

Different sets of wave functions are available in the literature. Here we use the functions 


derived in | 20 [] and adopted in JHJ to compute the coupling Qb*b-k- At the 6 -quark mass scale 
the values of the various functions appearing in eq.( [2.15| ), at u 0 = |, are: ^(wo) = 1-219, 
(p a (u. o) = 1.463, gi(u 0 ) = —3.37 x 1CT 3 GeV 2 and G 2 (u 0 ) = —1.8 x 1CT 2 GeV 2 ; these two last 
numbers correspond to the choice 5 2 {rrib) = 0.17 GeV 2 . The expressions for the functions 0 


and 0 can be found in | 1 (J; their contribution to the sum rule corresponds to a small fraction 
(of the order of 10 ~ 2 ) of the final result. 


Using mb = 4.6 GeV [^2|, < uu >= —(240 MeV ) 3 and the threshold So in the range 
s 0 = 33 — 35 GeV 2 we get the result depicted in Fig.2. The duality domain in M 2 , where 
the continuum contribution is small (less than 30 %), can be chosen as M 2 = 6 — 12 GeV 2 . 
The result is: 


/l* 9b*b*tt = (8.5 ± 1.5) x 10 “ 3 GeV 2 , (3.1) 

where the uncertainty is due to the variation with the threshold so and the Borel parameter 
M 2 . We do not include the uncertainty on the wave functions, which we consider as external 
input parameters. 

Eq.(|2.15|) can also be used to compute the coupling gD*D*rr- Using m c = 1.35 GeV 
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and So = 6 — 7 GeV 2 , taking into account the variation of the wave functions with the 
renormalization scale, we derive the result in Fig.3. In the range M 2 = 2 — 4 GeV 2 we get: 

/£, 9d* D*n = (1.8 ± 0.3) x 10“ 2 GeV 2 . (3.2) 

The results (ED. (13) favourably compare with the outcome of short-distance QCD sum 
rules f 2 B »g B * B ** = (9.4 ± 1.8) x 10~ 3 GeV 2 , f 2 D «g D *D** = (1.7 ± 0.4) x 10~ 2 GeV 2 . 

The analogous couplings B*Bn, D*Dtc , computed using the same set of parameters, are: 


f B f B • SB-B, = (7.9 ± 1.0) x 1(T 3 GeV 2 
/d/d* 9d*d.t — (1-5 — 0 . 2 ) x 10 2 GeV 2 


(3.3) 

(3.4) 


to be compared with the short-distance based calculation [1]]]: f B f B * 9b*Btt = (7.0 ± 1.5) x 
10 " 3 GeV 2 and fo.fo* 9 d*d* = (1.2 ± 0.3) x 10 “ 2 GeV 2 . 

The analysis of the asymptotic sum rule ( |2.16|) allows to derive g and, by taking the 
logarithmic derivative in 1/E , the parameter A. Using x/q = 1.2 — 1.4 GeV and E = 
0.8 — 2.0 GeV, we get, from the curves depicted in Fig.4: 


A = 0.40 ± 0.07 GeV 
F 2 g = 0.031 ± 0.006 GeV 3 ; 


(3.5) 

(3.6) 


the analogous result obtained in [|TUJ is F 2 g = 0.035 ± 0.008 GeV 3 . 

The calculation of the couplings can be done, using eqs. (|3.1| - |3Tl ), once the leptonic con¬ 
stants are known. Since the correlator ( f 2 . 1 |) has been computed at the leading order in a Sj it 
is consistent to use the leptonic constants determined at the same order. Then, one can use 
23|: f B * = 190 ± 10 MeV , f D * = 220 ± 24 MeV , f B = 150 ± 20 MeV , f D = 170 ± 10 MeV , 


with the results: 


9b*b*tt — 0.24 ± 0.05 
g B * Bn = 0.28 ± 0.05 

To derive the value of g , in the limit mj, 


9d*d*it — 0.31 i 0.08 

(3.7) 

9d*D n = 0.40 ±0.07 . 

(3.8) 


oo, one can linearly extrapolate the above results: 
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(3.9) 


gp*p*7T — g (l H-) 

m Q 

„ b \ 

gp*Pi: — g (l H-) , 

m Q 


(3.10) 


obtaining: 


g = 0.21 ±0.07 
g = 0.23 ±0.08 


a = 0.6 ±0.9 GeV 
b = 1.0 ± 1.0 GeV 


(3.11) 

(3.12) 


from gp*p* n and gp*pir, respectively. The results (|3.11|) , ( |3.12|) show that the same value 
for g is obtained by extrapolating to m b —> oo the values obtained for gp*p* n and gp*pn- 
One can conservatively quote g = 0.22 ± 0.10 as the result of the extrapolation. It is also 
worth observing that, although the fit does not accurately fix the parameters a, b, the 1 /rriQ 
corrections are rather sizable in the case of the charm quark. 

Let us now consider the calculation of g from eq. ( |3.6| ). Using the determination in [|24 
F = 0.30 ± 0.05, obtained neglecting radiative 0(a s ) corrections, one gets: 


g = 0.34 ±0.10 , 


(3.13) 


a result which, although higher in the central value, is compatible with the outcome of the 
fit. 


It is interesting to consider the inclusion of 0(a s ) corrections in the lcptonic constants, 
to investigate at least partially the role of the perturbative corrections in the determination 
of g. Using, therefore, fp* = 213 ± 34 MeU, fu* = 258 ± 26 MeU, fp = 180 ± 30 MeV, 
fp = 195 ± 20 MeV [j23| , one gets: 


gB*B*7 r = 0.19 ± 0.04 , gn*D*ir = 0.23 ± 0.05 

gB*B it = 0.21 ± 0.06 , g D*Dn = 0.30 ± 0.06 , 


i.e.: 


g = 0.17 ±0.06 , a = 0.6 ± 0.8 GeV 
g = 0.17 ±0.08 , b = 1.0 ± 1.4 GeV 
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(3.14) 

(3.15) 


(3.16) 

(3.17) 









The inclusion of the 0(a s ) terms in the leptonic constants does not spoil the prediction 
that the same g is obtained from the B*B*tt and B* Bn channels. The numerical value of g 
results modified at the level of 30%. 


Several determinations of the coupling g can be found in the literature 10,9,15,16,25,26 


with values varying up to g — 1. The analysis reported here points towards small values of 
g , in agreement with the outcome of the relativistic potential model [^I|, which shows that 
the nonrelativistic result g = 1 is reduced (g ~ 1/3) if a relativistic treatment of the light 
quarks is adopted. 


IV. LIGHT-CONE SUM RULES FOR g Bl B 0 n AND g' 

The same method described in the previous sections can be used to determine the cou¬ 
pling g', which weightes in the lagrangian ( |1.1|) the strong interaction between positive parity 
heavy mesons and pions. For a finite mass of the heavy quark, the coupling g Bl B 0 ix can be 
obtained from the correlator of a scalar Jg = qQ and an axial J p = <77^75 Q current: 

G^(p,p + q) = i j dx < n~(q)\T[d(x)b(x)^u(0)]\0 > e ipx 

= G®(p + q) li + GWq li (4.1) 

since the invariant function gets contribution from the poles Bq and B\. Defining 

< ^(Q)Bo(p)\B 1 (p + q,e) >= V 77 ^ mB i g BlBon e - q (4.2) 

with m Bl and m Bo the masses of the positive 1 + and 0 + states, respectively, and 

< 0| bu\B 0 >= (4.3) 

m b 

< 01>= fsjnB^v, , (4.4) 

we get the Borel transformed sum rule for g Bl B 0 n'- 
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2/ Bl 4,/a,4 a = j _ c -..,- )m! ^, w 

/TT^fe I 


+ g-^/M 2 


M 2 


fnl^l 


(1 - - 6{ £™' md) K(“o) - (> - “») 


+ m rf 


du 


f m 2 

2m b f 7T (l - u 0 )g 2 (u 0 ) + — —— -2mlip a (u 0 ) 


6 (m u + m d )' 


1 


- ^4m^/^[c/i(M 0 ) + G 2 (u 0 )] 


(4.5) 


ruo rl—ai 

m b f n dai da 3 (>~.<p\\+<p\\+2(p±-2<p 1 )(a 1 ,l-a 1 -a3,a3) 

JO JuQ—OLl 


2 M 2 /s, 


lUQ—ai 

r ' 1_ “° da 3 


lo a 3 

r u 0 r-l-ax da 3 


(ps-iriuo, 1 — Wo — O' 3, 03) 


dan / —1 — «i — «3, 03) 

wo - C ^3 


In the infinite heavy quark mass limit, it is easy to derive from ( |4.5|) the sum rule for g', 
using the scaling law /b x = ^=: 

2 < uu > 


F +2 g’ = j (i - e V ^2Ef^ n (u 0 ) + ^ — ^ > ^a(w 0 ) - M^o) + G 2 («o)] 


3 ./, 


E 


(4.6) 


E is the Borel parameter; A' is the binding energy for the s 1 = | + doublet: A' = mgj 0 —m b , 


and it can be obtained by the logarithmic derivative of (fh6|). 

Before performing the numerical analysis of eq. ([4.5|) , we notice a potential numerical 
problem related to the presence of terms of opposite signs and similar sizes contributing 
to the sum rule. An example is represented by the first two terms in ( |4.5| ), respectively 
proportional to tp n and ipp. For intermediate values of the heavy quark masses there is no 
evidence of a hierarchical structure among the contributions from different twists; moreover, 
cancellations between various terms results in a critical dependence on the values of the 
wave functions, which have their own uncertainties. On the other hand, in the infinite 
heavy quark mass limit such a problem does not occur, and therefore eq. (|4.6| ) provides us 
with the possibility of determining g’ more rcliaby. 

Using rriB 1 = rriBo = 5.7 GeV and so = 36 — 38 GeU 2 we get 


/bi/bo 9 b 1 b 0 tt — (1-2 ± 0.4) x 10 3 GeV 


(4.7) 
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A determination of gp 1 p Q7V immediately follows, using the leptonic constants of B\ and B 0 
computed in [[23|]: f Bl = 180 ± 30 MeV, f Bo = 180 ± 30 MeV (the corresponding constants 
in the charm case are f Bl = 240 ± 20 MeV and fn 0 = 170 ± 20 MeV). We get: 


g Bl Bon = 0.05 ± 0.02 . (4.8) 

Going to smaller heavy quark masses, the coupling gp lBon decreases. However, the cancel¬ 
lation between different terms becomes more and more important, and the result becomes 
critically dependent on the values of the wave functions. 

The asymptotic sum rule (|4.6|) , depicted in Fig.5 and obtained varying y 0 in the range 
yo = 1.3 — 1.5 GeV, and the Borel parameter E in the range E = 1.0 — 2.0 GeV, provides 
the result: 


A' = 0.89 ± 0.05 GeV 
F +2 g' = 0.022 ± 0.003 GeV 3 . 

Using F + = 0.46 ± 0.06 GeV^ P3j the last equation corresponds to: 


(4.9) 

(4.10) 


g' = 0,10 ± 0.02 


(4.11) 


The comparison between (|4.8|) and ([4-111) gives an insight on the role of the heavy quark 
mass corrections in the case of the determination of the strong coupling g'. 


V. LIGHT-CONE QCD SUM RULE FOR g BlB *-K AND h 

To complete our analysis of the couplings appearing in (|1.1| ) we consider the parameter 
h corresponding to the couplings BiB*tt and B 0 Btc. The coupling g BoB7T has been computed 


by light-cone and short-distance QCD sum rules in |TIJ, while a determination of g BlB * n can 


be found in [|lj]]. Here we want to derive an expression for h. 
Let us define g BlB *n by the matrix element 


< B*(p,e)n(q)\B l (p +q,Vj) >= g BlB *„(e* 


h)— + (p ■ v){q ■ ; 

m B , 


(5.1) 
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gB 1 B*-K is related to the coupling h in eq. (|l.l| ) by: 


2 h 


gB\B*Tv — — ~rV m Bi m B* , 

I TV 


(5.2) 


whereas the structure F is subleading in the uiq —> oo limit. To calculate gBiB*n one has 
to consider the correlator of an axial and a vector current: 


H^(p,p + q) = i J dx < 7T (g)|T[d(a;)7^6(a;),6(0)7 I/ 75M(0)]|0 > e ipx 
= H^g^ + H^ PlxPv + ... 


(5.3) 


The 1 + and 1 poles contribute to the invariant function H^°\ which must be considered to 
derive gBiB*ir- One obtains the Borel transformed sum rule: 


“r* + m s. 


gB 1 B*n { m Bi — m B*)fB 1 m B 1 fB*'< r l'B*e 2M 2 — 


M 4 , , , rribrrii 


( v /v 2 m u + m d 

2 M 2 (l + — [G' 2 (m 0 ) + g[ (u 0 )] 


M\ P [u 0)) - e~ m V M * 


2mlg 2 ( u o) 


+ M 2 


m 

M 2 
r u 0 da 3 


(5.4) 


c u 0 


'-Jo 2a 3 

da 3 

aai 


{<P\\ + <P||)(m 0 - «3, 1 - M 0 , 03 ) + 


da? 


'0 2 a? 


(<P|| - ^||)K, 1 - Mo - «3, Ol 3 ) 


^||(«1,1 - Oil - oi 3l a 3 ) 


IU0 — Ot.\ OL 3 

Differently from the calculation of g and g', the determination of h does not concern a 
symmetric correlator (Fig.l). However, also in this case the choice Uq = 1/2 allows us to use 
the same values of the pion wave functions employed for g , g', which reduces the uncertainty 
related to the wave functions. Considering that = ^(l) = G" 2 (^) — g [(|) =0, in 

the limit nib —> oo one obtains from eq.( |5.4| ) the expression: 

flrni E 

J 7T 7T 


F + Fh = 


-(1 — e y °/ E )e ip P (u 0 ) 


m u + m d A' — A 

Eq.(|5.5|) coincides with the expression that can be obtained from the sum rule for gB 0 Bn 


(5.5) 


11 


The numerical results from eq. (|5.4|) , using the same set of parameters chosen in the 
previous Sections: 


fB 1 fB*gB 1 B*n — 1-9 ± 0.4 GeV 2 , fD 1 fD*gD 1 D*-tr — 1-2 ± 0.2 GeV 2 


(5.6) 


17 


















produce the following values for gp^*^: 


9b x b*-k = 56 ± 15 , g Dl D*n = 23 ± 5 . (5.7) 

The linear extrapolation to mb —■► oo: hp 1 p* n = h( 1 + ^-) determines the coupling h: 

h = —0.7 ± 0.3 , d = 0.0 ± 0.7 GeV . (5.8) 


The same result comes from eq.( |5.5D ; as a matter of fact, choosing y Q = 1.3 — 1.5 GeV and 
the mass parameters A, A 7 derived in the previous sections, the result (Fig.6): 


FF+h = -0.083 ± 0.005 GeV 6 


(5.9) 


corresponds to 


h = -0.60 ±0.13 , 


(5.10) 


in agreement with the determination h 


—0.56 ± 0.28 obtained in 


0 


from g Bo Bn- 


VI. CONCLUSIONS 

The strong couplings of heavy mesons to pions are expected to obey symmetry relations 
in the mg —> oo limit, when QCD displays invariance under rotations in the heavy quark 
spin and flavour space. In this limit, it can be shown that the number of independent 
couplings reduces to a set of effective parameters. 

We employed light-cone sum rules to evaluate the couplings gp*p*-,T, gp*p l7r , gp 1 p 0 -n-, 
(P = D, B) both for finite values of the heavy quark mass and in the infinite mass limit. 
This procedure allowed us to observe that the light-cone sum rule method reproduces the 
right mq scaling behaviour for the physical quantities. Moreover, we verified that the asymp¬ 
totic results agree with usual QCD sum rules based on the short distance expansion. This 
comparison enforces our confidence in the reliability of both the variants of the QCD sum 


18 





rule tecnique (e.g. short distance based and light cone sum rules), which are among the few 
available instruments to investigate non perturbative aspects of QCD. 
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FIGURE CAPTIONS 


Fig. 1 


The correlator in eq. (|2.1|) . 


Fig. 2 

Stability curve for /^* gs*B* n- The curves correspond to the threshold s 0 = 34 GeV 2 
(continuous line), s 0 = 35 GeV 2 (dashed line), s 0 = 36 GeV 2 (dotted line). 


Fig. 3 

Stability curve for /|,» gn*D*w The curves correspond to s 0 = 6 GeV 2 (continuous line), 
s 0 = 6.5 GeV 2 (dashed line), s 0 = 7 GeV 2 (dotted line). 


Fig. 4 

Stability curve for A and for F 2 g. The curves correspond to the threshold yo = 1.2 GeV 
(continuous line), y 0 = 1.3 GeV (dashed line), y 0 = 1.4 GeV (dotted line). 


Fig. 5 

Stability curve for A' and for F +2 g '. The curves correspond to y 0 = 1.3 GeV (continuous 
line), yo = 1.4 GeV (dashed line), yo = 1.5 GeV (dotted line). 


Fig. 6 

Stability curve for FF + h. The curves correspond to y 0 = 1.3 GeV (continuous line), y 0 = 
1.4 GeV (dashed line), y 0 = 1.5 GeV (dotted line). 
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